00 

o 

O 
(N 

Sh 
Oh 1 

< 

00 
(N 

S 1 

Oh. 



CZ3 

O 
>V 

Ph. 



> 

00 

o 

oo 

o 
> 

X 



Low-frequency plasma conductivity in the average-atom approximation 

M.Yu.Kuchie\E 

Department of Theoretical Physics, School of Physics, 

University of New South Wales, Sydney, 2052, Australia 

W.R.JohnsorQ 

Department of Physics, 225 Nieuwland Science Hall, 

University of Notre Dame, Notre Dame IN 46556 

(Dated: April 28, 2008) 

Low-frequency properties of a plasma are examined within the average-atom approximation, which 
presumes that scattering of a conducting electron on each atom takes place independently of other 
atoms. The relaxation time r distinguishes a high-frequency region ujt > 1, where the single-atom 
approximation is applicable explicitly, from extreme low frequencies lut < 1, where, naively, the 
single-atom approximation is invalid. A proposed generalization of the formalism, which takes into 
account many-atom collisions, is found to be accurate in all frequency regions, from u) = to lot > 1, 
reproducing the Ziman formula in the static limit, results based on the Kubo-Greenwood formula for 
high frequencies, and satisfying the conductivity sum-rule precisely. The correspondence between 
physical processes leading to the conventional Ohm's law and the infrared properties of QED is 
discussed. The suggested average-atom approach to frequency- dependent conductivity is illustrated 
by numerical calculations for the an aluminum plasma in the temperature range 2-10 eV. 

PACS numbers: 51.70.+I, 52.20.fs, 52.25.Mq 



I. INTRODUCTION 

Various theoretical approaches are available to investi- 
gate the frequency-dependent conductivity of plasmas, 
ranging from methods based on a many-body expan- 
sion of the grand canonical partition function [l|, y, [3j 
to methods based on molecular dynamics simulations 
[J, \Mi Is 0, H M- I n tne present paper, we re-examine 
an average-atom approach [lOl ] that has been used re- 
cently to investigate anomalous dispersion in C, Al, Ag, 
and other plasmas in the soft x-ray region (14-47 nm) of 
the spectrum [0, GJ, [H, Q3, Q3] . The utility of the aver- 
age atom method rests on its simplicity and wide range 
of applicability. 

At low frequencies, electron-ion scattering contribu- 
tions dominate the conductivity <x(oj), while at higher 
frequencies (e.g. in the x-ray region mentioned above) 
photoionization and bound-bound transitions provide the 
most important contributions. Effects of multiple scat- 
tering were omitted in evaluating free-free contributions 
to a{uj) in Ref. [10| . leading to a (spurious) second-order 
pole at lo = that was regularized in an ad hoc way. 
In the paragraphs below, we discuss the origin of this 
pole in more detail and a give a modified formula for the 
free-free contribution to cr(uj) that accounts for multiple 
scattering, is regular at u> = 0, and rigorously satisfies 
the conductivity sum rule. 

The present discussion concerns the plasma conduc- 
tivity o~{oj) at low frequencies, i.e. presuming that the 



frequency to is lower than both the plasma frequency and 
typical frequencies of atomic excitations 



(j < (47m e e 2 /™) 1/2 , me 4 /h 3 



(1.1) 
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The only parameter that drives the conductivity in this 
region is the relaxation time t, which establishes a 
boundary between relatively high-frequencies lot > 1, 
which will be called the high-frequencies for short, and ex- 
treme low-frequencies, where lot < 1 including the static 
limit oo = 0; these frequencies will be called the ultra-low 
frequencies. 

Physical processes, which govern the conductivity in 
these two regions, differ qualitatively, as discussed in de- 
tail below. Alongside this physical difference, there ex- 
ists also a distinction in theoretical methods. One line of 
research is based on the Ziman formula, which is applica- 
ble in the static limit, leading to the conventional static 
Ohm's law, see Ziman [l6| ch.7, Mahan [17[ ch.8. An 
alternative approach is based on the Kubo-Greenwood 
formula [18|, [l9|, [20, [2l|, [22| , which usually gives a reliable 
description of conductivities at high-frequencies. 

Generally speaking, the Kubo-Greenwood formalism 
should lead to accurate results for arbitrary frequencies, 
provided though that all important scattering processes 
are taken into account. However, typically, within some 
given theoretical scheme, it is feasible to account only for 
some particular class of scattering events. This restric- 
tion may substantially reduce an area of applicability for 
the Kubo-Greenwood formalism. In particular, it is usu- 
ally difficult to extend its validity to the static limit. An 
example illustrating the latter fact is provided by models 
based on the average-atom approximation. In such mod- 
els, scattering of conducting electrons is assumed to take 
place on each atom, independently of other atoms. This 



means that many-atom events (multiple scattering), in 
which several atoms produce a coherent contribution are 
neglected. The simplicity and clear physical nature of 
the average-atom approximation make it popular. Its 
origins can be traced to the Thomas-Fermi model of 
plasma devised more than a half century ago by Feyn- 
man, Metropolis, and Teller [23[. A quantum mechanical 
version of the average-atom model is given by Blenski 
and Ishikawa [2J] and a recent implementation is found 
inRef. [H. 

The present work shows that scattering processes 
which take place at ultra-low frequencies necessarily in- 
clude several atoms. We will call these processes, the 
many- atom collisions for short. Their importance indi- 
cates that for ultra-low frequencies the single-atom ap- 
proximation breaks down. This fact explains a difficulty 
that occurs in the Kubo-Greenwood formalism in the 
static limit. The breakdown of the single-atom approx- 
imation manifests itself as a divergence of the conduc- 
tivity calculated in the Kubo-Greenwood formalism in 
the limit to — ► 0. As mentioned in the introduction, the 
conductivity in this limit exhibits a second-order pole. 

Thus, direct numerical calculations based on an 
average-atom approximation and relying on the Kubo- 
Greenwood approach are applicable for high frequencies 
only, while for lower frequencies, where the many-atom 
collisions are important, the formalism faces a difficulty 
This work resolves this difficulty, proposing a new ap- 
proach that is applicable for frequencies that satisfy the 
inequality (|1.1|) . In the static limit ui — 0, our descrip- 
tion reproduces the Ziman formula. For high frequen- 
cies ojt > 1, our results agree with the conventional 
Kubo-Greenwood description. In the intermediate re- 
gion lot ~ 1, the validity of our formalism is supported 
by the fact that it provides the correct result for the con- 
ductivity sum-rule. One of the important advantages of 
the proposed description is related to its simplicity. We 
show that all necessary physical quantities can be eval- 
uated using a simple single-atom approximation. This 
means that multiple scattering, which is paramount in 
the static limit, is accounted for effectively in the single- 
atom approximation! 

There is an important relation between the divergence 
in the conductivity at ultra-low frequencies and the in- 
frared problem of QED. To make this point more trans- 
parent, let us keep in mind that the conductivity de- 
scribes absorption and emission of quanta of the electro- 
magnetic field, which are possible due to electron scatter- 
ing. Presuming that the potential, which is responsible 
for scattering, is localized in a vicinity of some atom, one 
can express the amplitude of absorption / a b s (or emis- 
sion) in terms of the elastic scattering amplitude /. This 
relation reveals that the absorption amplitude has a pole 
at lo = 



f f 

Jabs OC — , 







(1.2) 



a photon line inserted into the outer electron legs as is 
shown in Fig. [1] The first-order pole in the absorption 
amplitude in Eq. (|1.2p leads to a second-order pole in the 
conductivity, 

<j(w) oc — . (1.3) 

LO A 

Developing this argument, we will show below that many- 
atom collisions prevent the divergence of the scattering 
amplitude in Eq. (|1.2[) at w — > 0. This happens because 
many-atom collisions lead to a finite relaxation time r, 
which measures the interval of time during which the 
electron travels between two subsequent collisions with 
different atoms. We show that the relaxation time pro- 
vides an effective cutoff for the amplitude in Eq. (ll.2|) . in 
which the pole is replaced by a finite quantity |/ a bs| — 
l/w m in = t. The well defined, finite, scattering ampli- 
tude leads to a conductivity that is regular at lo = 0; 
Eq. (|1.3p is replaced by the relation (j{lo) oc l/w^ in = t 2 . 
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This general, well-known, feature of the infrared pro- 
cesses in QED, is described by Feynman diagrams with 



FIG. 1: Two Feynman diagrams represent the amplitude, 
which describes electron scattering with absorption of a quan- 
tum of the electromagnetic field, which has frequency lo and 
small wave vector that leaves the electron momentum un- 
changed. The solid lines show the electron propagation, the 
dashed line - the quantum of the electromagnetic field, the 
filled dot - the elastic scattering process. In diagrams (a) and 
(b) the lines that represent the electromagnetic quantum are 
inserted into the outer legs, which makes these diagrams in- 
frared singular, oc 1/lo when lo — > 0. Other possible diagrams 
have no such singularity. 

Our discussion below presumes that the plasma con- 
ductivity is due mainly to scattering of conducting elec- 
trons by atomic cores. There exist other mechanisms 
contributing to conductivity. One of them is related to 
electron-electron scattering. The main idea of this work 
can be generalized to cover this mechanism (and others) 
as well. However, in order to keep our presentation sim- 
ple and clear, we will not attempt to formulate the idea 
in the most general case, restricting our discussion to 
electron-atom scattering only. At sufficiently high tem- 
peratures, when atomic cores are highly ionized, one ex- 
pects that electron- atom scattering gives the dominant 
contribution to the conductivity owing to the fact that 
scattering by an ion is a coherent process, with probabil- 
ity proportional to Z? on , where Z- lon is the ionic charge. 
By contrast, electron-electron scattering is an incoher- 
ent process with probability proportional to Z- mn . Thus, 
scattering by an ions is expected to dominate electron- 
electron scattering, provided Z\ QVl > 1. 



II. SINGLE-ATOM APPROXIMATION FOR 
HIGH FREQUENCIES 

A. Absorption of photons and elastic scattering 

We need to recall several simple important facts related 
to absorption of low-frequency quanta by electrons. Let 
us presume that there is a localized potential U = U(r), 
which causes electron scattering. Let us assume further 
that there is some external low-frequency homogeneous 
electric field. [25[ Then scattering can be accompanied 
by absorption of a quantum of the electromagnetic field. 
The process of absorption is described by the matrix el- 
ement / abs , 



/abs = (ipf\e-p\ipi) 



(2.1) 



Here e and p are the polarization vector of the electro- 
magnetic quantum and the operator of momentum of the 
electron, tpi and tpf are the wave functions of the electron 
in the initial and final states. We are interested in the 
low frequency region specified by Eq. ljl.ip . Our first goal 
is to simplify the matrix element in Eq. (|2.1|) , presenting 
it as 



considered as a diagrammatic representation of Eq. (|2.4|) . 
The singular energy denominator 1/oj, which appears in 
Eq. (|2.4p , arises only in the external legs of the two Feyn- 
man diagrams shown in this picture. All energy denomi- 
nators of all other diagrams include virtual energies of the 
atomic excitations, which are sufficiently high compared 
with the energy of the electromagnetic quantum. Corre- 
spondingly all other diagrams, which are not shown in 
Fig. [U are all finite in the limit u> — > 0. This fact distin- 
guishes the two diagrams in Fig. [TJ and guarantees that 
Eq. (|2.4p is accurate for low frequencies. 
From Eq. (|2.1[) one finds 



mLU 



(2.5) 



Averaging over possible orientations of the polarization 
vector one writes 

{(e-qf) = \ q 2 = l i {\- C os9) P \ (2.6) 

where brackets ( ) refer to the averaging procedure and 
9 is the scattering angle. Eqs. (|2.5l2.6[) give 



/abs = ^ (Vy I u\ p) - ^-(p'I u\ il> p ) (2.2) 

Here the first and second terms correspond to the Feyn- 
man diagrams, in which the line representing the elec- 
tromagnetic quantum is inserted into the left and right 
legs of the diagram respectively, see Fig. [JJ The wave- 
functions |p) and \ifip) in Eq. (|2.2p describe the electron 
propagation in the plane wave approximation and with 
account of the potential U, respectively. All processes, 
in which the line representing the electromagnetic quan- 
tum is inserted into internal parts of the diagram have no 
poles in the limit u — > 0, allowing one to neglect them in 
Eq. (|2.2|) (see the more detailed discussion after Eq. (|2.4p V 
Generally speaking, the electron energy in the initial 
and final states of the photo-absorption process are dif- 
ferent. However, for low frequencies this difference is 
insignificant. Neglecting it, one can presume that the 
matrix elements in Eq. (|2.2p are related to elastic events, 



i.e. |p| = |p'|. Remember now that the elastic scattering 
amplitude is defined as 

/=-^2^p'I C/ Ip) = -^2(p'I [/ IV'p), (2-3) 
27m 2-kTi 



(2.4) 



Consequently one finds from Eq. (|2.2p that 



Jabs 



2irh 



(e-q)/ 



where q = p' — p is the transferred momentum. 

Eq. (|2.4p relates the amplitude of the process with ab- 
sorption of a low- frequency electromagnetic quantum and 
the amplitude of elastic scattering [compare Eq. (|1.2D ]. 
Relations of this type provide a basis for the known in- 
frared problem in QED, see e.g. [2g|. Fig. [JJ can be 



J(\U s \ 2 )dn = l 



2 (2-Khp 



(l-cos8)\f\ 2 dn. (2.7) 



where the integration runs over the angles fi of the scat- 
tered electron. The factor 



<rtr= I {l- cos e)\f\ 2 <m 



(2.8) 



represents the transport cross section on the potential U. 
Eqs. (|2.7l2.8p give 



(\Us\ 2 )dn 



2 (2ixTiv 



a tT 



(2.9) 



Here v — p/m is the velocity of the electron. The quan- 
tity on the left-hand side of Eq. (|2.9|1 describes the proba- 
bility of absorption of low-frequency quanta. The trans- 
port cross section on the right-hand side is related to elas- 
tic scattering. A close connection between low-frequency 
electromagnetic processes and elastic scattering is well 
known, see e.g. Ref. [26j |. 

Deriving Eq. (|2.9p . we assumed that the potential U 
responsible for the electron scattering is localized within 
some finite volume. Precisely this property allows one to 
distinguish the two Feynman diagrams in Fig. [TJ Other- 
wise, if the potential is spread all over an infinite volume, 
a mere concept of an external leg of the diagram would 
make no sense. 

We can specify the potential U assuming that it is 
created by a single atom. In that case Eq. (|2.9p describes 
those events that take place during electron scattering 
by a single atom, being thus closely related to the single- 
atom approximation. 



B. Kubo-Greenwood formalism 

Consider the conductivity of plasma, which is due tc 
scattering of conducting electrons by atoms. Within the 
Kubo-Greenwood formalism it can be written as 



a(u>) = 



27m., e 2 



<</y|e-«|Vp>l (fp-V) ( 2 - 10 : 
d 3 p d 3 p' 



xS(s p -e p i -u) 



{2nh) 3 (2nh) 3 



where n a is the density of atoms, and f p is the Ferm 
distribution function for conducting electrons (which wil 
be denoted by f below) 



fr, 



1 



exp[(e p - fj)/kT] + 1 



(2.11 



The chemical potential /j, here is related to the concen- 
tration of conducting electrons n c 



d 3 p 

{2nh) 3 



(2.12) 



where the coefficient 2 accounts for two projections of 
spin. It should be noted that we have omitted contri- 
butions to Eq. (|2.10p arising from atomic bound states. 
These contributions, which lead to bound-bound reso- 
nances and singularities near photoionization thresholds, 
are insignificant in the low-frequency region of concern 
herein. 

The first factor in the integrand in Eq. (|2.10j) can be 
conveniently rewritten with the help of Eqs. (|2.1l2.9|) ; the 
difference of the distribution functions in the integrand 
can be simplified using the low frequency approximation. 
These transformations allow one to simplify the expres- 
sion for conductivity Eq. (|2.10| l. reducing it to 



a{u) 



2 n„.e 2 



3 uj 2 



v 3 a tT 



df\ d 3 p 



de) (2irh) 3 



(2.13) 



It is convenient to introduce the relaxation time r p for 
conducting electrons, which is due to collisions with 
atoms 



1 



V 71 a cr t r 



(2.14) 



Clearly, it depends on the electron momentum via the 
velocity and the transport cross section. Eq. (|2.13| ) can 
be written in this notation in a transparent compact form 



a(w) 



2e2 
3 



<9f\ d 3 p 



de) (2-Khy 



(2.15) 



This is, in fact, the low-frequency limit of the Kubo- 
Greenwood formula in the single-atom approximation. 

There are two conditions that restrict a region of fre- 
quencies in which Eq. (|2.15p is valid. Firstly, as was men- 
tioned, the frequency must be sufficiently low. More pre- 
cisely, this condition implies that the relevant scatter- 
ing phases 5i (s) , where I is a typical orbital momentum, 
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FIG. 2: The conductivity of an aluminum plasma at T — 5 
eV. Black line - the Kubo-Greenwood formula Eq. (l2.10|l . 
Green line - simplified Kubo-Greenwood formula Eq. (12.15|l . 
which predicts a second order pole cj{lu) oc 1/lu 2 . Red line 
- results from Ref. [10( | . which used an interpolating proce- 
dure to extend the results of the Kubo-Greenwood approach 
to the static approximation described by the Ziman formula 
Eq. pTTj) . Blue line - prediction of Eq. (|4~g)l . 



should not reveal significant variation in the interval of 
frequencies u> 



hio 



de 



<1 



(2.16) 



Here e is a typical energy of conducting electrons, 
ondly, the frequency is restricted from below by 

LJT V > 1 , 



Sec- 



(2.17) 



p being a typical momentum of those electrons that give 
significant contributions to the conductivity in Eq. (|2.15[) . 
We will discuss this condition in detail after Eq. (l4.8|) . 
Here, let us mention briefly that the necessary high fre- 
quency specified by Eq. (|2.17p makes it certain that scat- 
tering processes on different atoms take place incoher- 
ently, as independent events; in other words, that the 
single-atom approximation is valid. Eq. (|2.15[) predicts a 
simple oc 1/uj 2 behavior of the conductivity on frequency. 
If one ignores the restriction given in Eq. (|2. 1T[) by taking 
the static limit in Eq. (|2.15p naively, then this equation 
clearly indicates that the conductivity has a second-order 
pole at u = 0, as seen in Eq. (|1.3p . Fig. [2] illustrates this 
statement by comparing calculations based on the com- 
plete Kubo-Greenwood formula Eq. (|2.10p with predic- 
tions of Eq. (|2.15p . As an example, an aluminum plasma 
at temperature of 5 eV was taken in Fig. [5] The agree- 
ment between the two sets of calculations, shown in the 
black and green lines in Fig. [21 supports the validity of 
the approximations, which led to Eq. (|2.15p . The numer- 
ical code for calculations reported in the present work 



are based on the average-atom model and the numerical 
methods suggested in Rcf. (Iq |. 



III. CONDUCTIVITY AT ULTRA-LOW 
FREQUENCIES 

A. Ziman formula 

Consider the static limit lo = 0. The Ziman formula, 
which describes the conductivity due to electron-atom 
scattering reads 



a(0) = 



2e 2 



V 2 T„ 



~d~e 



d 3 p 

{2nh) 3 



(3.1) 



One observes its drastic distinctions from the result of the 
Kubo-Greenwood - type approach. Firstly, the Ziman 
formula Eq. (|3.1|) gives a constant static limit for the con- 
ductivity, while the Kubo-Greenwood formula Eq. (|2.15|) 
diverges at w = 0. Secondly, in these two formulas the 
conductivity shows an opposite dependence on the re- 
laxation time; Eq. (|2.15p reveals an inverse dependence, 
oc l/r p , while the Ziman formula (|3.ip predicts a linear 
dependence, oc r p . 

To find an origin for these distinctions, let us note that 
deriving Eq. (|2.15p we assumed that the electron momen- 
tum p is a good quantum number, which is changed only 
due to scattering on one given atom. Generally speak- 
ing, this assumption is incorrect. The momentum can be 
changed due to scattering on other atoms as well. 

To see the implications of this fact more clearly, let us 
note that the second-order pole oc 1/uJ 2 in the conductiv- 
ity arises as a direct consequence of the first-order pole 
oc 1/lu in the amplitude. The latter can be written as an 
integral 



1 _ _. 

LU 



,o 

/ exp(— itut) dt 

J — oo 



(3.2) 



where \t\ gives a period of time, which precedes an elec- 
tron collision with the given atom. Eq. (|3.2[) shows that 
deriving Eq. (|2.15p one presumes that during all this pe- 
riod of time, which can be very large, up to infinity, the 
electron momentum remains constant. 

As a matter of fact, this is not true. The momentum 
can remain constant only over a finite period of time, 
which equals a typical interval of time between two sub- 
sequent collisions. This interval is measured by the relax- 
ation time. This means that the relaxation time should 
necessarily produce the cut-off for the integral over time 
in Eq. 



/0 rO 

cxp(-iujt) dt — » — i / cxp(— iLut) . (3.3) 
-OO J —Tp 

The cut-off procedure can be fulfilled slightly differently 
and more conveniently, by introducing the cut-off func- 



tion in the integrand 

pQ p0 

— i\ cxp(—iojt)dt^—i / cxp(— iuit — |£|/r p ) 



1 



W + iJT p 



(3.4) 



The pole at uj = 0, which exists on the left-hand side 
here, is replaced by a finite behavior of the right-hand 
side. 

These simple arguments show that the pole oc 1/lu 
in the scattering amplitude and, correspondingly the 
second-order pole oc 1/lu 2 in the conductivity, are closely 
related to the single-atom approximation. The many- 
atomic events lead to the relaxation time, which erases 
this pole behavior. This argument is developed below, in 
Section HVl 



IV. 



RESONANT STATES OF CONDUCTING 
ELECTRONS 



According to Sec. IIII A[ multiple scattering events 
should play an important role in the ultra-low frequency 
region. In order to account for this scattering, let us start 
from a simple physical picture. If the electron has the mo- 
mentum p, then it keeps this momentum only for some 
finite period of time (relaxation time) because collisions 
with atoms in a plasma inevitably change it. In the classi- 
cal approximation this implies that only some finite part 
of the classical trajectory of the electron can be described 
by the initial momentum, while longer parts of the trajec- 
tory "forget" this momentum. Similarly, in quantum de- 
scription the stationary quantum states, which describe 
the electron propagation in a plasma, cannot be charac- 
terized by the momentum. 

However, if the plasma is sufficiently transparent, i.e. 
the relaxation time is sufficiently large, then during long 
intervals of time the electron momentum on classical tra- 
jectories remains constant. Consequently, the quantum 
states, which describe the electron propagation in a re- 
gion outside atomic cores during moments of time sepa- 
rating consequent collisions, should look similar to con- 
ventional plane waves. The fact that collisions, which 
destroy the electron momentum, are essential can be ac- 
counted for by stating that a quantum state with the 
given momentum p exists only during a finite period of 
time. In other words, the electron wave function of the 
conducting electron outside the atomic core of some atom 
resembles a conventional plane wave, but with the restric- 
tion that it exists only during a finite period of time that 
equals the relaxation time. Presuming that the relax- 
ation time is large, one can say that this wave function is 
a quasistationary state, which is similar to a plane wave, 
but possesses a finite width T p defined by the relaxation 
time 



r r 



(4.1) 



This identity, combined with Eq. (|2.I4p . states simply 
that T p — 2hvn a <7t r , which makes sense. 

The arguments just presented show that the electron 
wave function outside the atomic core of some atom can 
be written in a form 



\p,t) =exp \i(p-r 



£ P t) - Y l 



(4.2) 



This simple wave function possesses important physical 
properties. Firstly, it is close to a plane wave. Secondly, 
its finite width accounts for multiple-scattering events; 
i.e. collisions with different atoms. The width of this res- 
onant state is described by the relaxation time Eq. (|4.f[) ; 
the larger is the relaxation time, the closer is the wave 
function to a plane wave; exactly what one should expect 
when collisions are rare. 

The above argument can be developed further. If one 
wishes to consider the electron wave function in a close 
vicinity of a given atom, then the plane wave should 
be replaced by the wave function which takes into ac- 
count the influence of the atomic potential U(r). In 
other words, one needs to make in Eq. (|4.2[) a substitu- 
tion exp(ip ■ r) — > ip p (r). As a result, the wave function 
of a conducting electron, which takes into account the 
potential of a given atom, as well as scattering by other 
atomic particles, has the following form 



* p (r, t) = ip p (r) exp ( is p t - -^-t 



(4.3) 



Let us repeat, ipp( r ) here is the wave function, which 
describes the electron behavior in the potential created 
by a single atom, while T p is the width, which describes 
the momentum relaxation due to scattering processes on 
all atoms. 

It is instructive to compare Eq. (|4.3|) with a simple, 
classical idea of relaxation of the momentum. Consider 
for this purpose a value of the momentum averaged over 
the wave function Eq. 



(4.4) 



P(t) = ~J **(r,t)p*(r,t)d 3 



Here p is the operator of the electron momentum and V is 
a large, but finite volume, which makes a ratio in Eq. (|4.4|) 
well defined, V- independent. From Eqs. (|4.3| ). lj4.4|) one 
immediately finds that P(t) = exp(-T p t/K) P(0). The 
above can be written in a more routine form 



at t p 



(4.5) 



Clearly, this expresses a relaxation of the electron mo- 
mentum in conventional classical terms; t p plays here a 
role of the classical relaxation time, as one should have 
expected. Thus, a quantum description of the relaxation 
of the electron momentum based on the wave function 
Eq. (|4.3|) reproduces a well-known conventional classical 
physical picture. 



Using the wave function Eq. (|4.3[) in the Kubo- 
Greenwood formalism, one can follow a path outlined 
in Section III Bl However, a well-known short-cut makes 
these calculations redundant. A quasistationary nature 
of the wave function Eq. (|4.3p indicates that an ampli- 
tude of any resonant process involving this state acquires 
a conventional resonant energy denominator 



1 



AE + iTp/2 



(4.6) 



Here AE is a deviation of energy from its resonant value, 
which is presumed to be low. In our case this devia- 
tion is defined by the frequency of the electromagnetic 
field AE — fiuj. The resonant factor Eq. (|4.6[) coincides 
with the one found in Eq. (|3.4p . which underlines again 
a main physical idea; the multiple-scattering events al- 
low the electron momentum to exist only during a finite 
period of time. 

The resonant amplitude Eq. (|4.6[) always brings into the 
probability the resonant factor 



1 



AE 2 + T|/4 



(4.7) 



which is often called the Breit-Wigner factor (in atomic 
physics this describes Lorentzian lines). 

Applying Eq. (|4.7p to the process at hand one takes 



Eq. (|2.15p and, making the substitution l/o 
Tp/4), arrives to the following result 
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As mentioned earlier, the above result can be obtained 
directly using the wave function Eq. (|4.3p in the Kubo- 
Greenwood formalism; however, the abbreviated deriva- 
tion based on Eq. (|4.7p makes the discussion more trans- 
parent. 

Eq. (l4.8p differs from Eq. (|2.I5p . which was derived in 
the single-atom approximation, by the only physical fact; 
it accounts for many- atom collisions. This distinction 
becomes crucial in the static limit, allowing Eq. (|4.8p to 
reproduce correctly the Ziman formula Eq. (|3.f p . Thus, 
we return to the statement, which was mentioned sev- 
eral times previously; the many-atom events are very 
important for low frequencies. In contrast, in the high- 
frequency region ojt p > 1, Eq. (|4.8p reproduces Eq. (|2.I5p . 
which is based on the simple single-atom approximation. 
Thus, for ojTp > I the many-atom events become irrele- 
vant, in agreement with discussion in Sees. HI and III Bl 

Optical properties of a plasma are conveniently de- 
scribed with the help of a complex conductivity, which 
allows one to define the complex refraction index. Using 
Eq. (|4.8|) for the real part of the conductivity and apply- 
ing the conventional Kramers-Kronig dispersion relation 
Refs. [IllUmiHim, one finds that its real and imag- 
inary parts may be written side by side 
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Eq. (|4.8[) is one of the main results of this work. Its simple 
nature inspires a feeling that it could have been written 
without any discussion, as a simple convenient interpola- 
tion between the Ziman formula Eq. (|3.ip and the results 
of the Kubo-Greenwood approach Eq. (|2.15p . However, 
it is rewarding to realize that this result follows from a 
clear physical idea, which states that a conducting elec- 
tron can possess a constant momentum only over a finite 
period of time. 



V. SUM RULE 

Using Eq. (l4.8|) . one can calculate a simple but impor- 
tant integral 
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(5.1) 



Rewriting here d 3 p — m 2 v de p dQ and integrating over 
the energy by parts one finds 



a(u) = — e 2 m 2 / v 3 
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(5.2) 
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Taking into account Eq. (|2.12p , one finds that Eq. (|5.2|) 
represents the known, conventional conductivity sum- 
rule 



<j{w) dw 



n c e' 



(5.3) 



The integration in Eq. (|5.3p includes the region of very 
high frequencies, above the limit in Eq. (|l.l[) . which can- 
not be reliably covered by Eq. (|4.8[) . However, this region 
gives only a small contribution to the sum-rule because 
at high frequencies the integration in Eq. (|5.3|) converges 
very well, as J dcu/uj 2 . The sum-rule Eq. (|5.3p supports 
validity of Eq.gj]). 



VI. SUMMARY 

Let us summarize the main physical ideas. In a vicin- 
ity of a given atom the wave function of a conducting 
electron is strongly influenced by the potential of this 
atom. This fact makes it natural to presume that the 
problem can be formulated with the help of an average- 
atom model, which accounts for this variation and de- 
scribes correctly the electron scattering on this atom. 
However, the conductivity is related to processes of ab- 
sorption and emission of electric field quanta during scat- 
tering of the conducting electron by a given atom. When 
the frequency of the field is ultra- low, the absorption and 



emission take place in a region located far away from the 
atom. The electron wave function in this region is nec- 
essarily influenced by potentials of other atoms. As a 
result, the lower the frequency is, the more important 
are the many-atom events. Thus, the naive single-atom 
approximation inevitably breaks down in the static limit, 
where multiple scattering become crucial. 

From the first glance the necessity to account for mul- 
tiple scattering should make things much more compli- 
cated for the theory. There is though an important sim- 
plification. The many-atom events manifest themselves 
mainly via a restriction, which they put on the period 
of time during which the conducting electron can pos- 
sess a given momentum. Henceforth, one can account 
for these events by stating that the wave function of a 
conducting electron is a quasistationary state, which ex- 
ists only during a large, but finite period of time, which 
equals the relaxation time for the momentum. This idea 
can be expressed in terms of the quasistationary state 
Eq. (|4.3p . which describes the conducting electron. As a 
result, it becomes possible to account for multiple scat- 
tering within the formalism of the average-atom approx- 
imation, which greatly simplifies the problem. 

Applying this idea within the framework of the Kubo- 
Greenwood formalism, we find that the conductivity is 
given by Eq. (|4.8|) . which possesses several interesting 
properties. 

1. In the static limit it reproduces the Ziman formula 
Eq.(3l]). 



2. In the high-frequency region it is reduced to the 
Kubo-Greenwood type formula Eq. (|2.15p derived 
within the naive average-atom approximation. 

3. It satisfies the conventional sum-rule Eq. (|5.3p . 

4. It is formulated in terms of physical quantities, 
which can be evaluated in the average-atom ap- 
proximation that is convenient for applications. 



Our starting point was Eq. (|2.4p . which relates the elastic 
scattering amplitude with the amplitude of absorption 
(emission) of low frequency quanta. The latter gives a 
particular example of a general property of QED, which 
allows one to express any radiation process with soft 
quanta via a purely elastic scattering process, see e.g. 
[26j . Starting from Eq. (|2.4|) . we derive Eq. (|2.15p using 
a single-atom approximation, then taking into account 
many-atom events upgrade it to Eq. (|4.8p which, in the 
static limit, reproduces the Ziman formula Eq. (|3.ip for 
Ohm's law. Thus, the well known, conventional Ohm's 
law may be considered as a direct consequence of general, 
fundamental infra-red properties of QED. 

Using Eq. (|4.8p and evaluation the necessary average- 
atom quantities r p and f p with the help of the model of 
Ref. [10( we calculated the conductivity of the aluminum 
plasma at different temperatures. The results are shown 
in Figs. [21 and [3] In the one-atom approximation the 
conductivity is divergent as a oc 1/lo 2 , as shown in the 
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rely on an interpolation. It is satisfying that the two 
approaches give close numerical results; compare the red 
and blue lines in Fig. O 

Fig. [3] presents results of a series of calculations based 
on Eq. (|4.8p for different temperatures of the plasma. 
The increase of the conductivity with temperature re- 
flects an increase of the concentration of conducting elec- 
trons. The rapid decrease of the conductivity in the high- 
frequency region underlines the important role played by 
low frequencies; as was mentioned, the low-frequency 
region gives a dominant contribution to the sum-rule 
Eq.flOD. 

In conclusion, it is shown that Eq. (|4.8[) successfully 
describes plasma conductivity at low-frequencies. 



FIG. 3: The conductivity of an aluminum plasma at differ- 
ent temperatures. Calculations are based on Eq. (|4.8[l . The 
parameters used in this formula are evaluated with the help 
of the approach of Ref. |10l |. 



black and green lines in Fig. [21 To avoid this unphysical 
divergence Ref. [T(| suggested a particular interpolating 
procedure, shown by a blue line in Fig. [2J which brings 
the conductivity to the Ziman formulae Eq. (|3.1[) in the 
static limit. Eq. (|4.8p provides more rigorous treatment 
of the conductivity at low frequencies, which does not 
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